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NOTATION 


a5j^  “ 

- 

C  (Wave) 
C  2 (Wave) 


c  (SubBorlpt) 

D 

D  (Subscript) 
dyn  (Subscript) 
B  » 

"■sc' 

"'Tp” 

L  - 


Slant  distance  that  the  wave  propagates  (See  Figure  3*^) 

Slant  distance  that  the  wave  propagates  (See  Figure  3*^) 

Dllatatlonal  (Wave) 

Shear  (Wave) 

Velocity  of  wave 

Velocity  of  wave 

Velocity  of  sound  In  air 

...  of  the  shell 
diameter  of  the  shell 
...of  the  medium  (soil) 

Dynamic  value  of  . . . 

Modulus  of  elasticity 
Thickness  of  shell 
"Ambient"  constant 

Length  of  the  cylindrical  shell 


I  - 

m  (Subscript) 

P  - 

^0 

p  (Subscript) 
t  « 


Structxiral  deflection 
...In  the  model  (shock  tube) 
Overpressure 
Ambient  air  pressure 

. . .In  the  prototype 
any  time 
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MOTAnOW  (Continued) 


"A 

W 


t 

A 

t 


*82'*S1  * 


Time  when  C  wave  leaves  point  A  (See  Figure  3*4) 

X 

Time  when  C  wave  leaves  point  B  (See  Figure  3«^) 

2 

Time  when  C  wave  arrives  at  S  (See  Figure  3*^) 
a  1 

Time  i^n  C  wave  arrives  at  S  (See  Figure  3>^) 

2  2 

Pulse  duration 

Time  li^se  between  arrival  at  the  structure  of  C  and  C  waves 

i  2 

Time  lapse  between  the  two  parts  of  a  wave  split  by  a  non-linear 
material 


U  - 
u  “ 

V  *  u  » 


^4yn 


c 

-  -i 


Air  shock  velocity 
Displacement  of  the  medium  (soil) 
Particle  velocity 
Any  distance  coordinate 


Depth  of  burial 


dyn 


1  -  2i/, 


D,  d>n 


Height  density 
Flexural  rigidity 


Loearlttanlc  decreaunt  for  percyclic  dancing 
Strain 

Mach  angle  of  C  wave 
1 

Mach  angle  of  C  wave 

2 
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X 


A  = 
“ 

V  - 

P  “ 
a  - 

T  - 


WOTATlCTi  (Continued) 

Scale  factor  of  the  quantity  Indicated  by  its  subscript 
(See  Table  4.1) 

Wave  length 

Characteristic  length  of  response  wave 

Poisson's  Ratio 
Hass  density 
Stress 

Any  natural  period  of  the  structure 


Rote:  The  notation  appeeurlng  in  the  appendices  is  defined  where  used. 


1.0  innQiiicriOK 


Tbe  stndj  pxograa  pursued  imder  the  present  contract  has  been  directed 
tofiard  tbe  objective  of  initiating  a  logically  directed  approach  to  the  prediction 
of  soil  and  structure  response  to  air-blast  effects  of  ttaemonuclear  weapons  by 
■eans  of  experiaents  with  snail  nodels  subjected  to  slnulated  blast  loading. 
Specifically,  it  la  intended  to  ^  ^  tbeoreticall^r  based  foundation  for  the 

future  direction  of  aboek-tid>e  requlrenents. 

It  is  inportant  to  recognise  the  philosophy  that  nust  tmderlle  the  approach 
to  the  preaent  stikhr.  Generally,  a  vide  gap  exists  between  theory  and  practice 
in  tbe  use  of  nodels  to  stu^y  physical  phenosKna.  Ibis  is  nalnly  due  to  tbe 
difference  in  interests  of  the  theoretician  and  the  ea^rinenter.  Because  of 
this,  nodel  eaqierlaents  have  not  always  been  as  fruitful  as  could  otherwise  be 
eaqpected.  The  prdhlens  of  nodel  ea^rinentatlon  cannot  be  successfully  solved 
exclusively  by  theory  or  experinent,  but  1y  a  caqAtlble  cohbination  of  both. 

liblle  there  has  been  effort  eaq^ended  in  both  theory  and  eii^rlnent  in 
previous  nodel  studies  in  the  area  of  underground  structure  response  to  alrblast 
loadings  no  unifying  approach  has  been  advanced.  The  present  study  is  directed 
toward  gathering  the  Infomatlon  produced  in  previous  nodel  studies  and  in 
theoretical  work  and  welding  it  into  a  framework  upon  which  future  model  exper- 
Ixents  my  be  based.  Such  future  e]q;ieriaents  then  will  be  both  directed  and 
purposeful. 

Thus,  the  present  study  is  the  first  link  in  a  balanced  overall  progran  of 
theory  and  e:q;ierinent.  The  results  of  this  first  stage  will  not  be  a  panacea 
for  nodel  test  procedures,  but  will  be  tbe  basis  for  tbe  initial  i^iase  of 
meaningful  eq^rlxentatlon. 
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The  philosophy  of  the  present  study  then  is  one  of  perspective.  It  gathers 
experience  and  theory  from  the  past  and  applies  It  to  the  production  of  more 
meaningful  results  In  the  future.  It  will  mainly  determine  the  questions  to  he 
answered  hy  future  work  and  the  proposed  methods  for  obtaining  these  answers. 

1.1  Objectives 

The  fundamental  problem  with  which  this  study  program  Is  concerned  Is  the 
development  of  scaling  laws  which  will  enable  the  stinictural  response  of  under¬ 
ground  structures  subjected  to  air-blast  loading  to  be  predicted  by  means  of 
model  experiments. 

This  problem  Involves  the  application  of  dimensional  analysis  to  formulate 
a  means  for  qualitative  (and,  where  possible,  quantitative)  prediction  of  full- 
scale  behavior  of  dynamically  loaded  structures  through  interpretation  of  model 
experiments. 

The  problem  also  Involves  the  use  of  distorted  scaling  and  model  simulants 
where  such  are  necessary  to  achieve  the  desired  end  result.  Determination  of 
the  type  and  usefulness  of  the  Information  that  may  be  gathered  from  model 
experiments  is  another  facet  of  the  problem.  The  results  of  the  study  program 
consist  of  recommended  methods,  procedures,  and  scaling  laws  for  future  shock- 
tube  experiments. 

1.2  Approach 

The  research  program  under  the  present  contract  has  been  divided  into  three 
phases;  (l)  Preliminary  Studies,  (2)  Development  of  Scaling  Law,  and  (3) 

Model  Material  Investigations. 

Preliminary  studies  Include  a  literature  search  and  ascertainment  of  the 
characteristics  of  the  Air  Force  shock  tube. 
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During  the  conduct  of  the  study  progran  the  approach  to  the  second  phase 
vas  naterially  altered.  Initially,  as  reported  in  the  Interla  rqport,  this  phase 
vas  subdivided  into  a  well  ordered  stystem  of  categories  of  distinct  phenomena. 

These  phenosiena  (loading,  soil,  structure,  action  of  loading  on  soil,  and  inter¬ 
action  between  soil  and  structure)  and  their  attending  parameters  were  first 
considered  more  or  lees  iu  vacuo.  Then  it  vas  jdanned  to  use  the  relationships 
developed  for  these  distinct  catagories  and  the  classical  methods  of  dimensional 
analysis  and  theory  of  models  to  formulate  scaling  relationships.  The  conflicts 
encountered  between  the  scaling  relationships  thus  devdoped  were  to  be  adjusted 
and  resolved  into  a  meanlngfful  scaling  law. 

As  this  itese  of  the  worh  continued  it  became  apparent  that  the  approach 
was  leading  to  a  coagdexly  intertwined  and  conflicting  scaling  ^stem. 

Tundamental  pAienoaKna  and  peripheral  trivia  were  so  inextricably  coaibined  that 
the  basic  questions  to  be  answered  ly  eiqperlmental  processes  were  obscured. 

At  this  Juncture,  a  re-evaluation  of  the  problem  led  to  a  revision  of  the 
method  of  approach.  It  was  decided  to  concentrate  on  the  aforementioned  basic 
questions  and  to  determine  how  they  could  best  be  answered  in  as  unfettered  a  manner 
as  possible. 

In  accordance  with  the  revised  approach,  maximum  effort  vas  directed  toward 
determining  which  questions  conceiving  the  general  problem  of  dynamically  loaded 
underground  structiures  were  most  in  need  of  answers.  Extant  theories  were  re¬ 
viewed  to  determine  upon  what  assunq^tlons  their  major  premises  were  based.  Con¬ 
flicting  theories  were  coaqiared  to  determine  wherein  th^  differed  and  how  they 
could  be  resolved.  These  considerations  led  to  the  formulation  of  several  pro¬ 
positions  which  were  used  as  the  basis  for  the  scaliig  relationships  developed. 
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This  approach,  vtaen  taqilaDented  hy  shock  tube  eaperineotatlon,  vUl  serve 
to  either  conf Im  or  disprove  current  leading  theories  in  the  fld.d  of  dynamically 
loaded  underground  structures.  Such  confliaation  or  denial  should  lead  to  a 
■ore  rapid  advance  in  theory  and  thus  to  a  new  set  of  shock  tube  eaqieziinents. 

This  then  is  in  conforaance  with  the  phllosqpfay  of  this  study  as  presented  in 
Section  1.0, 

Phase  3  of  the  prograa  consisted  of  developing  an  outline  of  the  aethod  to 
he  used  for  determining  the  properties  of  aodel  siaulants  necesscuy  to  properly 
scale  the  desired  effects.  A  list  of  possible  simulants  and  their  properties 
was  also  prepared. 

In  the  derivation  of  the  scaling  laws,  the  shape  of  the  structure  used  was 
a  cylinder  in  all  cases.  This  allowed  the  derivation  of  the  scaling  law  to  be 
as  uncluttered  as  possible  and  will  make  e3q;>erimentation  easier.  (VHoders  are 
relatively  easy  to  fabricate  (as  opposed  to  spheres,  domes,  and  the  like).  Since 
the  phenomena  of  Interest  are  general  rather  than  directed  to  given  structures, 
a  single  geometiy  will  serve  the  purpose  of  indicating  the  effects  of  the  phen¬ 
omena  Instead  of  the  peculiarities  of  different  shapes. 
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2.0  SIATEHERF  OF  THE  FRCfflLEK 


In  order  to  retain  perspective,  the  statement  of  the  problem  vlU  he  In 
three  parts  corresponding  with  the  three  phases  of  the  program  as  outlined  In 
Section  1.2.  This  Is  not  to  la^y  that  the  aforementioned  3  phases  are  co-equal 
In  tu^ortance  or  that  th^  represent  equal  amounts  of  effort.  The  Important 
phase  and  that  upon  which  the  most  effort  was  expended  was  the  second.  Hie 
great  hulk  of  this  report  will  deal  with  phase  2. 

2.1  Preliminary  Studies 

Ttie  first  step  In  conducting  an  Investigation  of  the  use  of  models  to 
simulate  4)naa>dcally  loaded  imderground  structures  was  to  determine  the  results 
of  previous  Investigations  pertinent  to  the  problem.  To  accomplish  this,  the 
problem  was  divided  Into  three  sections  and  a  literature  search  was  conducted  to 
establish  the  present  state  of  the  art  In  each  field.  The  areas  selected  for  the 
literature  search  were: 

a.  Blast  slBif  list  Ion  techniques  (shock  tube) 
h.  IVnaadc  behavior  of  soils 

c.  Theoretical  Investigations  involving  the  response  of  burled  shell- 
type  structures  to  a  shock  wave. 

The  literature  survey  accaqilished  its  intended  purpose  of  Indicating 
present  approaches  to  the  problem.  A  bibliography  of  the  portion  of  the  liter¬ 
ature  surveyed  which  was  most  helpful  to  the  present  problem  Is  presented  In 
the  Appendix  of  this  report. 

A  visit  to  the  Air  Force  Shock  Tube  Facility  and  discussions  with  the  operat¬ 
ing  personnel  thereof  served  to  acquaint  AMP  personnel  with  the  physical 
realities  of  the  model  experiment  progpram  which  will  follow  the  present  study. 
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2.2  Development  of  Scaling  Lav 


In  order  to  develop  a  scaling  law  that  was  not  just  another  rehash  of  the 
familiar  applications  of  the  theory  of  models  as  developed  for  problems  in 
hydraulics  and  aerodynamics,  it  was  necessary  to  formulate  a  fresh  approach  to 
the  specific  problem  at  hand-how  to  make  use  of  models  to  determine  the  pheno¬ 
menology  of  dynamically  loaded  underground  structures.  Unfortunately,  the  use 
of  models  was  conceived  in  the  field  of  hydraulics  and  developed  in  the  field  of 
aerodynamics.  This  genesis  has  led  to  stereotyped  applications  of  the  theory 
of  models  to  other  problem  areas  and  provided  a  road  block  to  those  endeavoring 
to  use  models  in  other  fields.  Thus,  it  was  decided  to  forego  tradition  in  the 
present  study  in  order  to  develop  a  scaling  law  which  would  be  meaningful  to 
the  solution  of  Important  problems  in  the  field  being  investigated. 

Essentially,  the  procedure  Involved  the  formulation  of  basic  conditions 
that  reflected  the  important  problems.  Then  the  scaling  law  was  developed. 

2.2.1  Formulation  of  Basic  Conditions 

The  experience  of  MRD  personnel  in  the  field  of  underground  structures  and 
blast  loading  and  the  results  of  the  literature  search  have  led  to  the  formulation 
of  a  list  of  basic  conditions  or  questions  which  bear  heavily  on  the  xiltimate 
solution  of  the  problem  of  buried  structures.  The  entries  on  this  list  are  not 
necessarily  co-equal  in  significance  or  importance. 

Two  objectives  were  used  to  guide  the  formulation  of  this  list.  The  first 
was  to  determine  the  relevance  of  certain  phenomena  whose  significance  to  the 
problem  was  of  a  debatable  nature.  The  second  was  to  develop  the  capability  of 
obtaining  at  least  qualitative  information  as  to  certain  phenom.'n  .Ich  are  felt 
to  be  lii5)ortant  but  have  not  yet  yielded  to  theoretical  Investigations, 
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The  list,  as  formulated,  heginB  with  the  importance  of  the  role  of  structural 
response  in  the  total  problem.  Under  what  circumstances  does  the  response  of  the 
structure,  when  combined  with  the  nature  of  the  forces  exerted  upon  the  structure, 
lead  to  one  history  of  deformation  and  under  what  set  of  circumstances  does  it 
lead  to  another?  The  first  family  of  conditions  deals  with  the  relation  of  the 
natural  period  (or  periods)  of  the  structure  and  the  pulse  width  of  the  incoming 
wave.  In  postulating  this  family,  the  ratio  of  the  pulse  width  to  a  natural 
period  of  the  structure  was  investigated.  The  range  of  values  which  this  ratio 
could  run  was  divided  into  3  regimes.  Regime  I  includes  those  values  of  the 
ratio  where  the  pulse  width  is  extremely  short  in  comparison  with  the  natural 
period  of  the  structure.  Regime  II  contains  pulses  whose  durations  are  roughly 
one-half  as  large  as  that  of  the  natural  period  of  the  structure. 

Regime  III  Includes  the  remaining  higher  values  of  the  ratio.  The  boundaries  of 
these  regimes  have  been  worked  out  and  this  analysis  is  included  in  later  portions 
of  this  report. 

The  second  condition  developed  in  the  consideration  of  a  realistic  scaling 
law  also  involved  the  natural  periods  of  the  structure.  This  condition  arose 
from  consideration  of  the  question:  Which  modes  of  the  structure  that  are  ex¬ 
cited  by  the  dllatatlonal  wave  are  of  such  frequency  and  mode  shape  that  they 
are  reinforced  by  the  shear  wave  by  virtue  of  its  time  of  arrival?  If  such 
relationships  do  exist,  it  was  felt  that  they  must  be  preserved. 

Consideration  of  the  influence  of  Mach  angles  led  to  the  third  condition 
Imposed  on  the  scaling  system.  If  the  Mach  angles  must  be  preserved,  the  ratio 
of  the  shock  velocity  in  the  air  to  the  seismic  velocity  must  be  preserved  since 
it  is  this  ratio  which  determines  the  Mach  angle.  However,  it  was  determined 
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(a) 


FLane  Haves  Striking  a  cyiin<3er 
Tigure  2.1 


8 


that  the  requirement  to  preserve  Mach  angles  is  not  a  strict  one.  Consider, 
for  example,  a  structure  of  the  form  of  a  horizontal  cylinder  whose  longltud- 
irial  axis  is  parallel  with  the  wave  front.  This  example  is  Illustrated  in 
Figure  2.1.  Geometrically,  since  the  cylinder  is  symmetric  about  its  axis,  the 
result  is  basically  the  same  in  Figures  2.1(a)  &  2.1(b).  The  limiting  case  is 
9  =  90“,  for  beyond  this,  the  superselsmlc  case  is  not  preserved.  Thus,  on 
purely  geometric  considerations,  M'-ch  angles  do  not  have  to  be  precisely  pre¬ 
served.  The  only  consideration  would  be  to  preserve  superseismicity.  It 
should  be  pointed  out,  however,  that  the  degree  of  non- preservation  of  the  Mach 
angle,  even  within  the  broad  limitation  imposed,  does  have  an  effect  on  the 
ratio  of  stress  Inqplnglng  on  the  structure  to  overpressure  at  the  surface.  As 
the  Mach  angle  Increases,  this  ratio  decreases.  These  considerations  then  form 
another  condition  for  the  scaling  law  developed. 

It  is  Interesting  to  note  here  the  relative  insensitivity  of  the  Mach  angle 
to  changes  in  air- shock  velocity  under  certain  conditions.  This  may  be  seen  in 
Figure  3*5  which  shows  Macb  angle  as  a  function  of  shock  velocity  for  different 
values  of  seismic  velocity.  Consider,  for  example,  the  curve  for  c  =  I5OO  ft/sec. 
Between  values  of  10,  CXDO  and  20,000  ft/sec  for  shock  velocity,  the  Mach  angle 
varies  only  from  about  9*  to  4®.  Such  behavior  can  be  useful  in  the  practical 
solution  of  difficult  scaling  problems. 

Non-linearities  of  the  soil  stress- strain  curve  formed  the  basis  of  the  next 
condition.  Presuming  that  soil  exhibits  a  stress-strain  curve  like  that  shown 
in  Figure  2.2,  the  questions  that  arise  concern  the  need  for  preserving  the 
shape  of  this  curve  in  the  model. 
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Figure  2,2  Typical  Stress-Strain  Curve  for  Soil 


For  example,  assuming  that  the  scaling  law  developed  would  be  non-linear  and 
that  a  and  e  would  be  scaled  Independently  of  each  other,  the  question  of  the 
relationship  between  the  loading  and  unloading  portions  of  the  curve  could  be 
significant,  since  the  slope  of  the  unloading  portion  of  the  curve  affects 
recovery  displacement  and  thus  energy  propagation  away  from  the  structure. 

Another  condition  imposed  was  that  involving  the  relationship  between 
soil  displacement  and  structural  deflection.  Since  soil-structure  Interaction 
depends  in  part  on  the  relative  motions  of  both  soil  and  structure  and  the  main¬ 
tenance  or  non-maintenance  of  contact  between  the  two,  it  becomes  necessary  to 
scale  soil  displacement  and  structural  deflection  the  same. 

In  addition  to  the  above  major  basic  conditions  the  following  factors  also 
influence  the  derivation  of  a  scaling  law. 

Depth  of  burial  of  the  structure  affects  the  problem  in  several  ways.  It 
may,  for  example,  affect  the  amount  of  attenuation  of  the  stress  wave  if  the 
attenuation  is  due  to  plastic  or  percyclic  dissipation.  It  is  also  involved 
in  determining  whether,  in  the  case  of  shallow  burial,  apy  significance  may  be 
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assigned  to  the  reflection  of  stress  waves  from  the  ground  surface.  Depth  may 
also  be  significant  if,  again  in  the  case  of  shallow  burial,  the  mass  of  soil 
above  the  structure  is  significant. 

When  a  structure  has  a  high  acoustic  impedance  compared  with  that  of  the 
soil,  the  effect  of  the  response  of  the  structure  may  become  significant. 

(structure  responses  can  cause  attenuation  of  local  soil  pressure  which  propa¬ 
gates  as  what  is  termed  the  "response  wave.")  If  these  phenomena  are  significant 
(this  is  probably  dependent  on  the  soil  recovery  capacity)  the  characteristic 
length  of  the  response  wave  must  be  scaled  as  the  incoming  wave  length. 

The  analysis  of  the  above  conditions  and  other  limiting  factors  will  be 
found  in  Section  3* 

2.2.2  Derivation  of  Scaling  Relationships 

The  end  product  of  the  aforementioned  analyses  was  expressed  in  terms  of 
scaling  relationships.  Certain  scale  factors  were  deemed  fixed  or  given  (such 
as  overpressure  and  therefore  shock  velocity)  and  the  remaining  scale  factors 
were  expressed  in  terms  of  these.  Each  desired  condition  dictated  a  number 
of  such  relationships.  These  were  tabulated  in  a  form  such  that  the  scaling 
law  for  any  particular  combination  of  desired  conditions  may  be  readily  ascertained 

From  the  scaling  relationships  and  associated  tables,  the  phenomena  which 
formed  the  bases  for  the  desired  conditions  can  be  incorporated  into  an  exper¬ 
imental  program  and  these  phenomena  can  be  imposed  on  model  structures.  By 
varying  the  values  of  the  scale  factors,  the  effect  of  the  phenomena  may  be 
studied  both  separately  and  in  concert.  In  this  way,  more  qualitative  (and  quan¬ 
titative)  information  may  be  obtained  to  enable  further  theoretical  studies  to  be 
made.  The  existence  of  other  effects  than  those  described  herein  may  also  be 
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tested.  ^  verifying  the  existence  or  non-existence  of  certain  phenomena,  a  clear¬ 
er  picture  of  the  total  problem  can  be  obtained. 

Suggestions  as  to  the  Implementation  of  the  above  approach  are  contained  in 
Section  6.0.  The  suggested  tests  given  there  are  by  no  means  all  inclusive,  but 
are  Intended  to  point  up  basic  problems  which  might  be  solved  via  the  experimental 
approach.  Other  problems  which  may  be  of  present  interest  to  AFSWC  personnel 
may  be  eunenable  to  this  approach.  It  is  also  expected  that  the  results  of  the 
testing  program  will  uncover  new  problem  areas  which  in  turn  can  be  made  the  sub¬ 
ject  of  new  tests. 

2.3  Model  Investigation 

After  the  scaling  laws  have  been  developed,  the  successful  application  of 
these  laws  will  depend  on  the  ability  to  find  model  materials  whose  physical 
properties  are  compatible  with  the  development  of  usable  scale  factors.  The 
material  parameters  of  a  large  group  of  available  and  inexpensive  materials  must 
be  obtained  either  through  published  data  or,  when  the  required  data  is  unavailable 
from  this  source,  by  tests  conducted  to  furnish  these  parsimeters. 

As  an  aid  in  choosing  model  material  simulants  several  important  criteria 
are  outlined  and  illustrated  by  an  example.  In  addition  a  tabulated  list  is 
Included  which  presents  some  possible  model  simulant  materials  and  their  material 
properties  which  are  pertinent  to  the  development  of  scaling  factors.  While  not 
all  inclusive  either  in  extent  or  pertinent  parameters,  it  is  Intended  to  furnish 
a  guide  for  the  initial  experimental  program  which,  it  is  expected,  will  follow 
the  theoretical  work  presented  in  this  report. 
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3.0  iwMtYsis  or  oCTromoais 


In  order  to  develop  a  set  of  scaling  relationships  based  on  the  afore- 
Mentioned  approach  of  setting  given  requirements  which  are  directed  toward  inm- 
ediate  problems,  it  is  necessary  to  analyse  the  conditions  postulated.  This 
section  is  concerned  with  the  analyses  of  these  conditions. 

The  coodltlons  are  divided  into  two  groups,  those  which  must  hold  for  those 
scaling  relationships  idilch  are  intended  to  explore  a  given  iriienomenon  and  those 
which  must  hold  for  all  the  phenomena  considered.  These  two  groups  will  be  con¬ 
sidered  separately.  Definitional  statements  which  serve  to  li^dement  the  analysis 
and  are  the  basis  for  developing  scaling  relationships  will  be  presented  where 
needed.  All  definitions,  definitional  statements>  ^aod  cohditlcihs  ubed  are  listed  in 
Section  k.  The  structure  chosen  for  all  cases  is  a  cylindrical  shell. 

3.1  Conditions  Imposed  on  All  Scaling  Relationships 

One  of  the  basic  prescribed  parameters  of  the  shock  tube  is  the  overpressure. 

TiV  means  of  the  Ranklne-Bugoniot  equation,  the  velocity  of  the  shock  front  may 


be  determined. 


-  {1  +  ff-) 


(3.1) 


given  as 


A  valid  approximation  of  this  equation  for  higher  overpressures  may  be 

z 

■  kp  where  k  (3.2) 


.  _  6  o 


Figure  3.1  shows  the  valid  range  of  the  approximation.  For  overpressures 
over  to  pel,  the  value  of  U  given  bX  equation  3.2  is  within  10](  or  less  of 
that  value  given  by  equation  3*1  >  lY  overpressiires  of  over  2^  psi  are  maintain¬ 
ed,  the  approximation  should  give  adequate  results. 

Because  of  the  boundary  condition  existing  between  the  air  and  soil,  the 
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vedue  of  stress  In  the  soil  at  the  surface  must  be  the  same  as  the  overpressure. 

For  this  reason,  the  condition  is  Imposed  that  the  scale  factor  for  stress  In  the 
soil  Is  the  same  as  the  scale  factor  fon  Overpressure*  (x  =  X  ) 

The  boundary  condition  between  the  soil  and  the  structure  Imposes  another  con¬ 
dition.  The  stress  in  the  soil  must  equal  the  pressure  on  the  structure  at  any  par¬ 
ticular  point.  This  does  not  mean  that  the  stress  field  in  the  simulated  soil  is  re¬ 
quired  to  be  the  scaled  stress  field  of  the  prototype. 

This  seme  boundary  condition  Imposes  the  condition  that  the  displacement  of  the 
structure  must,  for  any  given  point,  be  equal  to  the  soil  displacement  at  that  point, 
3.2  Conditions  Desired  for  Useful  Scaling 

The  different  phenomena  on  which  it  has  been  deemed  desirable  to  base  the  devel¬ 
opment  of  scaling  laws  are  hereinafter  analyzed.  Basic  formulae  upon  which  the  re¬ 
lationships  have  been  developed  are  contained  in  Tables  k,2  and  4.3* 

3.2.1  Pulse  Width  and  Natural  Period  of  the  Structure 

One  of  the  many  factors  which  determine  the  manner  in  which  a  structure  responds 
to  dynamic  loading  is  the  relation  between  the  pulse  width  of  the  loading  stress  and 
the  natural  period  of  the  structure.  Thus,  in  the  case  where  the  pulse  is  extremely 
short  in  comparison  with  a  natural  period  of  the  structure,  it  may  be  regarded  simply 
as  an  Intpulse.  On  the  other  hand,  in  the  case  where  the  pulse  is  long  and  has  a  long 
rise  time  compared  with  a  natural  period,  structural  response  is  affected  by  maximum 
pressure  behind  the  wave  front  rather  than  duration  of  the  losid. 

In  the  transitional  case,  both  the  peak  pressure  and  the  duration  time  affect 
structure  response.  Rise  time  also  has  an  Independent  effect  on  response.  For  the 
applications  which  are  of  future  Interest,  however,  very  short  (near  zero)  rise  times 
dominate  consideration.  Therefore,  variable  rise  time  has  not  been  considered  in 
this  study, 
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Consideration  of  a  simple  mass- spring  system  with  various  time- dependent 
forcing  functions  has  led  to  a  delineation  of  the  three  cases  mentioned  above. 
Figure  3*2  gives  graphs  of  the  results  of  this  computation.  From  these  graphs, 
it  has  been  concluded  that  if  the  total  pvilse  duration  <  0.1  of  the  natural 
period  under  consideration,  the  pulse  may  be  treated  as  an  impulse.  If  the  pulse 
is  over  0.8  of  the  natural  period,  peak,  pressure  is  the  dominating  factor.  Graph 
(b)  of  Figure  3-2  shows  the  effect  of  rise  time,  but,  as  previously  mentioned, 
this  has  not  been  taken  into  account  in  the  derivation  of  the  scaling  law.  The 
following  table  gives  the  limits  of  the  regimes  and  the  scaling  conditions  per¬ 
tinent  to  each. 


Regime 

I 

II 


III 


TABLE  3.1 

Ratio  of  Pulse  Width 
to  Natural  Period _ 

0  to  0.1 


0.1  to  0.8 


0.8  Or  more 


Scaling  Condition _ 

Scale  structural  deflection  in 
accordance  with  impulse  criteria 

Preserve  ratio  of  pulse  width  to 
natural  period,  scale  structure 
deflection  in  accordance  with  peak 
pressure  criteria 

Scale  structure  deflection  in 
accordance  with  peak  pressure 
criteria 


In  regime  I,  the  table  shows  that  structural  deflection  must  be  scaled  in 
accordance  with  Impulse  criteria.  Thus 


(3.3) 
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Figure 


Responses  to  equivalent  Dirac  inqjulse 


Regime  II  has  two  conditions  which  are  pertinent.  The  first  of  these  states 
that  the  ratio  of  pulse  width  to  natural  period  must  he  preserved  from  model  to 
prototype .  Thus 


=  X. 


W 


(3.4) 


The  second  qualification  is  that  structure  deflection  must  he  scaled  according  to 
peak  pressure.  This  dictates  that 


(3-5) 


In  the  third  regime,  deflection  of  the  structure  must  he  scaled  according  to 
peak  pressure  only.  Thus  3 

^  ’"da 

X,  =  (3.6) 


3.2.2  Arrival  Times  of  and  Cg  Waves 

Considering  the  fact  that  there  are  two  wave  fronts  being  propagated  at 
different  velocities  through  the  medium  and  that  both  these  fronts  will  impinge 
on  the  structure,  the  question  arises  as  to  what  influence  their  combined  effect 
will  have  on  the  total  problem.  More  specifically,  which  modes  of  vibration  of 
the  structure  that  are  excited  by  the  dllatatlonal  wave  are  of  such  frequency 
and  shape  that  they  are  reinforced  by  the  shear  wave  by  virtue  of  its  time  of 
arrival?  If  such  a  phenomenon  is  important,  it  should  be  preserved. 

The  velocity  of  the  dllatatlonal  wave  is  of  course  greater  than  that  of  the 
shear  wave,  and  it  produces  a  different  particle  motion  within  the  medium. 

Figure  3-3  Illustrates  the  subject  phenomenon.  For  purposes  of  'liflcation, 
the  waves  have  been  shown  as  plane  waves  with  no  loss  in  general!  ,y.  The  dir¬ 
ection  of  propagation  of  the  waves  in  each  case  is  perpendicular  to  the  front. 
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structure  Sia>Jected  to  Cj^  &  Cg  Ifeves 
Figure  3.3 


^9 


The  particle  motion  contributed  by  each  wave  is  shown  by  arrows  on  the  illustration. 
In' (a),  the  waves  are  approaching  the  structure.  A  short  time  later,  the  dilata- 
tional  wave  impinges  on  the  cylinder  (b).  The  waves  continue  and  the  structure 
deflects  (c).  If  a  vlbj’ational  mode  of  the  structure  causes  it  to  have  a  cross- 
sectional  shape  as  shown  in  (d)  at  the  same  time  as  the  shear  wave  contacts  it, 
the  result  will  be  a  reinforcement  of  the  vibrational  mode. 

For  the  purposes  of  this  analysis,  the  velocity  of  the  shock  front  in  air 
(U)  is  assumed  to  be  constant.  This  assumption  is  very  reasonable  for  present 
purposes  since  the  length  of  the  soil  bln  over  which  the  wave  passes  is  small. 

If,  for  example,  the  overpressure  in  the  shock  tube  is  100  psi,  the  shock  velo¬ 
city  will  be  roughly  2800  feet/second.  Thus,  the  wave  will  pass  over  the  soil 
bin  in  around  .003  seconds.  There  will  be  very  little  change  in  shock  front 
velocity  in  this  short  time.  The  assiimptlon  goes  a  long  way  to  simplify  the 
analysis,. 

Considering  now  a  point  in  the  medium,  the  problem  is  to  determine  the 

difference  in  arrival  time  between  the  two  waves  (See  Figure  The  point 

in  the  medlxom  has  been  given  two  designations  (S^  and  Sg)  to  distinguish  between 

the  arrived  times  of  the  two  waves.  When  considering  the  arrival  of  the  first 

wave,  the  point  will  be  referred  to  as  S  ;  when  considering  the  second  wave,  it 

will  be  called  G  . 

a 

From  the  figure,  we  see  that 
sin  (90*  -  Qi)  = 

and 

sin  (90*  -  62)  =  z/^^ 


20 


s  s 


1  2 


Arrival  Times  of  &  Cg  Waves 
Figure  3-^ 
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Since 


sin(90-  -  eo  = 


and 


sln(90°  -  6z) 


_  (U^  - 


U 


(U^  -  _  z  _ 

U  “  °x(ts^  -  t^) 


and 


(U^  -  C22)V2  _ 

u 


(3.7) 


(3.8) 


(it  may  be  seen  from  the  figure  that  AS,,  and  BS„  are  equal  to  c  (tr,  -  t. ) 

i  *=  1  D  ^  A 

and  (tg  -  tg)  respectively.) 


One  more  relationship  is  needed  to  establish  the 
(ts  -  tg  ).  The  geometry  of  the  figure  establishes 


required  time  interval 
the  necessary  equation  as 
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(3.9) 
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By  algebraic  manipulation  we  may  obtain 


2  1 


y.. 


% 


(3.10) 


which  Is  the  delay  time  we  are  seeking. 

From  this  equation  we  may  detezmlne  that,  for  a  given  stock  velocity  and  soil 
(or  soil  simulant)  the  time  between  arrival  of  the  waves  Is  purely  a  function 
of  the  depth  of  burial. 

A  slsqpllfylng  approximation  that  Is  useful  In  ranges  when  U  c  Is 


t  -  § 


By  using  the  definition 
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the  following  Is  obtained 


) 


A 


(3.12) 


.  /2(1-V 

*.\/ _ Sx 

°  1  V 


T 

dyn 
^^,dyn 


(3,13) 


3.2.3  Mach  Angle  Preseirvatlon 

'As  was  stated  In  Section  2.2.1  (and  Illustrated  In  Figure  2.1)  the  require¬ 
ment  to  preserve  Mach  angle  Is  not  necessarily  strict.  What  Is  necessary  Is  that 
the  Buperselsmlc  case  be  preserved. 

If  the  Mach  angles  were  to  be  precisely  preserved,  the  sines  of  these  angles 
would  have  to  remain  Invariant  between  the  model  and  the  prototype.  Since  the 
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sines  of  the  Mach  angles  6^  and  62  are  equal  to  the  velocity  ratios  Ci/Uand  C2/U 
respectively,  the  resulting  relation  between  shock  velocity  and  Mach  angle  is  as 
shown  in  Figure  3*5  for  selected  values  of  seismic  velocity.  As  long  as  these 
ratios  remain  <  1,  the  waves  will  retain  their  required  superseismic  character. 

More  precisely,  since  the  dilatational  wave  has  the  greater  velocity,  if  we 
restrict  the  ratio  of  it  to  the  shock  velocity  to  values  less  than  one,  the 
shear  wave  is  insured  to  be  superseismic.  Therefore,  our  scaling  law  may  be 
flexible  in  that  the  ratio  (ci/u)  in  the  model  may  vary  from  something  more  than 
zero  to  almost  1«  As  this  ratio  in  the  model  varies  from  that  in  the  prototype, 
however,  the  ratio  of  stress  at  a  given  point  in  the  medium  to  the  overpressure 
at  the  surface  also  varies. 

3.2.4  Non-Linear  Stress-Strain  Curve 

Perhaps  the  greatest  difficulty  in  the  theoretical  solutions  of  underground 
structures  subjected  to  nuclear  blast  effects  lies  in  the  non-linear  character 
of  the  soil- stress- strain  curve.  Figure  2.2  gives  a  typical  example  of  a  stress- 
strain  curve  for  soil. 

A  great  deal  of  useful  information  can  be  derived  from  using  a  soil  simulant 
which  will  produce  a  similar  curve.  The  scale  factor  for  strain  will  most  prob¬ 
ably  be  different  from  that  for  stress,  ^y  obtaining  a  close  similarity 
of  the  curves  foi  soil  and  simulant,  the  scale  factor  for  slope  of  the  stress- 
strain  curve  would  then  be  the  ratio  of  the  scale  factor  for  stress  to  the  scale 
factor  for  strain. 

Perhaps  the  most  important  aspect  of  the  non-linear  stress- strain  curve  in 
regard  to  scaling  is  its  Irreversible  character.  TL's  should  be  p  I’ved  in  the 
soil  simulant.  Such  a  property  can  probably  be  simulated  best  with  foamed  materials. 
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Since  the  slope  of  the  unloading  portion  of  the  stress- strain  curve  affects 
recovery  displacement,  and  hence  energy  propagation  away  from  the  structure,  it 
would  be  advantageous  to  be  able  to  scale  this  the  same  as  the  loading  portion. 
This  puts  another  qualification  on  the  soil  simulant.  It  also  requires  more 
knowledge  of  the  actual  soils  which  are  to  be  simulated.  However,  as  more  of 
these  conditions  are  met,  more  qualitative  information  can  be  obtained  which  will 
lead  to  greater  advances  in  theory. 

If  the  non-linear  stress- strain  curve  can  be  simulated  as  indicated  above, 
another  phenomenon  can  be  scaled.  Consider  the  case  where  the  shape  of  the 
nonlinearity  is  such  that  it  splits  the  pvilse  in  two.  The  analysis  of  these  two 
waves  and  the  difference  in  their  times  of  arrival  at  a  given  point  is  not  unlike 
the  case  of  the  time  difference  between  the  Cj.  and  C2  waves  as  given  in  Section 
3.2.2.  Here  we  may  say  that 


'B 


%  substituting  the  equivalents  of  c_  and  c.  we  see  that 

15  A 


X  xVz  x/^2 
^  ^D.dvn 


(3.15) 


3y  requiring  the  scaling  of  the  jjeriod  of  the  structure  to  conform  with  the 
scaling  of  t,  the  effect  of  such  a  sjxLit  pulse  can  be  determined. 

3.2.5  Structural  Deflection  and  Soil  Displacement 

The  stress  wave  in  passing  through  the  soil  medium  generates  a  localized 
displacement  of  the  soil  particles  along  the  path  of  the  wave  front.  When  the 
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wave  front  contacts  the  surface  of  the  buried  structure,  the  displaced  soil  ad¬ 
jacent  to  this  surface  exerts  a  force  on  the  structure  wall.  This  force  results 
in  a  displacement  (£)  of  the  structure  surface  the  magnitude  of  which  is  depend¬ 
ent  on  the  magnitude  of  the  exerting  force  and  on  the  specific  structural  para¬ 
meters  which  are  pertinent  to  the  phenomenon.  Since  the  displacement  of  the 
buried  structure  wall  must  be  compatible  with  the  corresponding  displacement  of 
the  soil  adjacent  to  the  surface  of  the  structure  at  any  given  point,  it  is  nece¬ 
ssary  that  the  displacements  be  scaled  the  same  in  both  media. 

Assuming  that  the  buried  structure  is  a  thin  walled  cylindrical  shell,  it 

(49) 

is  shown  in  the  theory  of  cylindrical  shells  that  the  inextensional  defor¬ 
mation  of  a  circular  cylindrical  shell  will  involve  an  expression  for  the 
deflection,^,  of  the  form 


Op  L  D^d-Ve") 

Ec  h3 


(3. .16) 


It  is  seen  from  the  above  expression  that  the  structural  parameters  per- 

cf 

tinent  to  the  structural  deflection  are  E  ,  L  ,  It— I  ,  and  v  .  Soil  displace- 


is  given  by 


u  = 


(3.17) 


so  that  for  sceuLing  purposes  the  relationship  governing  soil  displacement  is 


given  by 


X,.  X. 

Pd 


(3.iB) 


The  desired  condition  is  that  both  stresses  and  displacements  (and  hence 
strains)  in  the  structure  and  adjacent  soil  have  not  only  point  to  point  equal¬ 
ity,  but  that  the  scale  factors  for  these  remain  invariant  throioghout  the  model. 
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3.2.6  Attenuation 

Attenuation  of  stress  waves  through  dissipative  mechanisms  of  the  medium 
is  an  important  consideration  in  the  problem  of  buried  structures  subjected  to 
blast  effects.  One  of  the  jxjssible  mechanisms  by  which  such  attenuation  occurs 
is  through  percyclic  damping. 

Percyclic  damping  by  definition  obeys  the  law  such  that  each  sinusoidal 
conqponent  of  a  propagating  wave  is  damped  by  a  constant  reduction  factor  in 
prox>agatlng  a  distance  equal  to  its  own  wave  length.  Thus,  a  waveform  would 
change  with  time  according  to  the  equation 


f 

f(x,t)  =J  e  g(0 


(3.19) 


where 


f(x,o)  =  j  6(0®  is  the  initial  shape  of  the 


1  **  Ifx 

waveform  and  S(0  =  S  f(x,o)  e  dx  is  the  spectrum. 


Thus  assuming  percyclic  damping  as  the  attenuation  mechanism,  the  ratio 
of  the  initial  wave  amplitude  to  the  amplitude  at  some  distance  ct  is  given 


6ct 

"IT 


6ct 


by  e 


Since  e  is  a  transcendent&il  function,  the  argument  -^,9.^.  must  be  dlmen- 

A 

6ct 

slonless.  Thus  to  preserve  similitude  in  the  danqplng  ratios,  ■  must 

remain  invariant  between  model  euid  prototype.  In  addition,  the  wave  length 
of  the  propagating  wave  must  be  scaled  as  the  depth.  This  requirement  need  only 


28 


be  met  where  percycllc  daii5)ing  ie  to  be  scaled 


3.2.7  Response  Wave 

The  responses  of  a  structure  always  cause  an  attenuation  of  local  soil  pressure. 
When  a  structure  has  a  high  acoustic  Impedance  in  con^rison  with  that  of  the  eoil> 
this  attenuation  comes  gradually  and  only  after  the  sharp  rise  upon  Initial  reflec¬ 
tion.  This  attenuation  propagates  away  as  what  is  termed  the  response  wave.  Any 
characteristic  length  of  the  response  wave  is  proportional  to  the  ratio  of  the  stinic- 
ture  mass  per  area  in  the  median  of  the  shell  to  the  soil  mass  per  volume.  In  order 
to  preserve  the  effect  of  this  phenomenon  in  the  model,  it  will  be  necessary  to 
scale  the  response  wave  the  same  as  the  incoming  wave. 

The  relative  importance  of  preserving  this  ratio  is  dependent  on  the  recovery 
capacity  of  the  soil  in  question. 

3.2.8  Engulfment  Time 

The  time  it  takes  for  one  wave  length  of  the  stress  wave  to  pass  over  the  struc¬ 


ture  is  of  some  Importance.  When  the  ratio 


D 


CiT 


is  in  the  neighborhood  of  1 


(e.g.  1/3,  1/2,  1  1/2  etc.),  the  scale  condition  should  be 


(3.20) 


1 

If  the  ratio  in  the  prototype  is  much  different  than  1,  this  ratio  need  not  be 
scaled  closely. 

3.2.9  Scaling  of  Weight 

It  is  Inescapable  that  weight  density  must  be  scaled  the  same  as  mass  density, 
at  least  as  long  as  we  are  earthbound.  Thus, 


(3.21) 


(3.22) 
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Since  weight  is  a  force,  it  appears  at  first  that  it  should  be  scaled  as  other 
forces.  This  scaling  would  actually  relate  stresses,  in  the  form 


(3.23) 


(3.24) 


Bu  ,  scaling  of  weight  as  force  might  conflict  with  other  scaling  relations.  If 
such  conflict  arises,  there  are  several  possible  remedies. 

It  is  likely  that  in  many  cases  the  stresses  due  to  structure  weight,  and 
possibly  soil  weight  also,  will  be  only  a  fraction  of  the  stresses  due  to  blast 
overpressure.  In  such  cases  3-24,,  and  possibly  3.23,  can  be  ignored. 

If  weight  is  responsible  for  a  major  portion  of  total  stress,  an  attempt 
should  be  made  to  deal  only  with  overpressure- induced  stress,  which  would  be  the 
excess  above  weight- induced  stress.  This  amounts  simply  to  a  translation  of  the 
stress  scale,  since  weight  remains  constant,  and  should  not  Interfere  with  the 
dynamic  phenomena  being  studied.  If  this  can  be  done,  then  3.23  and  3.24  can  again 
be  ignored.  Difficulty  might  be  encountered  in  translating  the  stress  scale  of 
the  soil  due  to  nonlinearity  of  its  stress- strain  curve.  The  stress- strain  curves 
of  the  prototype  and  of  the  model  must  now  be  similar  after  translation.  But 
chances  are  50-50  that  this  should  present  no  more  difficulty  than  if  they  had 
to  be  similar  before  translation. 

As  a  last  resort,  to  account  for  soil  weight, X  should  be  replaced  by  X 

P  p  +  Z7d 

ip  all  primary  scaling  relations  except  the  one  relating  overpressure  to  air- shock 
velocity.  If  the  effect  of  structure  weight  cannot  be  neglected,  it  can  always 
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be  eliminated  by  a  translation  of  the  scale  of  structure  stress  because  the 
stress-strain  curve  of  the  structure  is  linear.  Again  both  3*23  and  3*24  can  be 
ignored. 
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4.0  DERIVED  SCALING  RELATIONSHIPS 


Table  4,5  gives  the  derived  scaling  relationships.  The  scale  factors  are  given 
in  terms  of  known  or  prescribed  scale  factors  for  the  given  desired  conditions  set 
forth  in  Section  3-2  and  listed  in  Table  4.4.  The  imposed  conditions  (given  in 
Section  B-l  and  listed  in  Table  4.3)  are  included  in  the  derivations  when  applic¬ 
able. 

Two  of  the  desired  conditions  have  not  been  Included  in  the  derived  scaling 
relationships  -  the  condition  for  preserving  superseismicity  and  the  scaling  of 
weight.  As  was  pointed  out  in  Section  3.2.3>  the  preservation  of  the  superseismic 
case  is  not  an  exceedingly  strict  requirement  and  it  may  be  checked  after  the  rest 
of  the  scaling  requirement  has  been  developed  (i.e.  check  to  see  that 
If  weight  is  required  to  be  scaled,  the  procedure  given  in  Section  3-2.9  should 
be  followed. 

In  the  column  for  each  desired  requirement  in  Table  4.5,  certain  parameters 
are  listed  as  given.  Most  often,  the  given  parameters  Include  the  physical  para¬ 
meters  of  the  structure.  These  are  either  easily  controllable  or  dictated  by  the 
physical  limitations  of  the  shock  tube.  Other  given  i)arameterB  must  be  ascertained 
experimentally  or  Just  chosen.  An  example  of  the  former  category  would  be  the 

determination  of  X  ,  X„  .  and  X  .  Since  these  refer  to  dynamic  pro- 

Cl-  P.dyn 

parties  of  the  soil  and  the  soil  simulant,  they  cannot  be  computed  from  published 
static  values.  Known  field  test  data  must  provide  the  values  for  the  soil  and 
"free-field"  experiments  in  the  shock  tube  must  provide  the  values  for  the  soil 
simiHant.  In  cases  where  the  scale  factors  for  properties  of  the  structure  are 
considered,  the  model  material  will  have  to  be  chosen  before  the  other  scale 
factors  are  determined.  This  may  involve  several  attempts  at  arriving  at  the 
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TABLE  4.1 


SCALE  FACTOR  lEEWTIFICATlON 


Scale  Factor 


tlyn 


^T) 

^*SynfP 


\da\ 

I'Sel  ,<iyn 


,dyn,m 

^dyn,p 


Paraineter 

overpressure 

stress  In  soil 

strain  in  soil 

deflection  in  structure 

displaceaent  in  soil 

dlajneter  of  structure 

depth  in  soil 

length  of  structure 

lynaaic  soil  modulus 

tine 
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TABLE  4.1  (Continued) 
SCALE  FACTOR  IDENTIFICATION 


Scale  Factor 


^Wjin 

"^wTp 


Parameter 

pulse  width  (in  units  of  time) 

weight  density  of  structure 

weight  density  of  soil 


^(/3-l),dyn 


ratio  involving  seismic  velocities 


U 

U 


m 

P 


bending  stress  in  structure 


velocity  of  dilatatlonal  wave  in  medium 


density  of  medluiu 


velocity  of  shock  front  in  air 


any  natural  period  of  structure 


thickness  of  structure  shell 


density  of  shell  material 
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TABLE  4.1  (Continued) 

SCALE  FACTOR  IDENTIFICATION 

Scale  Factor 

Parameter 

II 

’V 

"ambient"  constant  for  air 

II 

E 

flexural  rigidity  modulus  =  - 

l-v/ 

II 

ratio  of  Dla.  of  stinacture  to 
thickness  of  structure 

^,m 

characteristic  length  of  response  wave 

simulant  which  will  give  usable  scale  factors. 


Each  column  in  Table  4.5  gives  the  scaling  relationship  imposed  by  a 
desired  condition.  Columns  A(i)>  A(2),  and  A(3)  are  mutually  exclusive  in 
ascertaining  a  scaling  law  for  a  given  experiment.  Column  E  need  only  be  Included 
when  it  is  desired  to  scale  percyclic  damping  and  column  G  need  be  included  only 
when  tne  conditions  as  cited  In  Section  3-2.8  warrant  it. 

To  ascertain  the  scaling  relationship  for  a  given  test,  the  conditions  to 
be  met  must  be  decided  upon.  These  conditions  dictate  the  scaling  factors. 

Where  two  or  more  relationships  for  one  quantity  are  encountered  (one  dictated  by 
each  separate  condition  which  involves  that  quantity),  these  must  be  made  compat¬ 
ible.  For  example,  in  order  to  have  conditions  B  through  G  hold  with  condition 
A(1),  must  be  scaled  as  dictated  by  condition  A(i)  and  can  no  longer  be  chosen 

freely.  By  the  same  token,  the  requirement  for  scaling  I  dictates  that 
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TABLE  4.2 


DEFINITIONAL  STATEMENTS 


=  /3  ^ 

c  dyn 

2 


D,  dyn 


T  = 


D^  P,V. 


h  A  Vs 


A  = 


1-v. 


t  = 


^  (^-1) 


dyn 


=  z 


Va 


Pt 


■/  .Vs 

/  \ 

do 

do 

'B,dyn~ 

1  de  / 

Vs 

A ,  dyn 


thus  restricting  the  choice  of  the  scale  factors  involved. 

The  choice  of  A(l),  A(2),  or  A(3)  depends  on  the  relation  of  the  pulse  width 

to  the  natural  period  as  explained  in  Section  3*2. 1.  The  time  sceile  factor  in  A(l) 

t 

is  given  by  Xu  =  where  and  t^  represent  measured  time  intervals  of  the  load¬ 

ing  wave  in  shock  tube  and  prototype  respectively.  Figure  4.1  shows  an  example 
of  this. 
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TABLE  4.3 


PRIMARY  RELATIONSHIPS  IMPOSED  FOR  ALL  CASES 


a 


(point  wise) 


Prototype  Shock  Wave 
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TABU  k»5 

DERIVED  SCALIHQ  RBLATIOWSHIPS 


Deaired 


Condition 

A 

B 

c 

Seale 

Factor 

(1) 

(2) 

«j| 

• 

* 

fl 

• 

« 

• 

« 

« 

* 

X* 

x^  H 

J  x»  xi 

• 

1 

D 

• 

"h 

1 

8 

a 

• 

« 

• 

• 

\ 

* 

X- 

P  t 

V.4 

i 

X^  X 

Pd  ®i 

xi  xi 

^Po 

Hfl 

X 

a  i 

X“X*  X* 

D  P„  P 

z 

- 5^ 

^(A-l),dyn 

X.  xj  X^  X 
^  ^  Pd  ®D,dyn 

^4 

_ 

*  Denota  "given**  or  **cho8en"  ecaie  factor  for  that  case 
I  0  &  h  are  here  given  only  by  ratio  O/h 

Blank  spaces  indicate  no  requirement  for  this  factor  In  the  particular  case 


For  all  cases  s  are  given  or  chosen 


X  .  x‘  x‘ 
\  *lc  P 


-  K 


4o 


The  pulse  width  (t  )  of  A(2)  may  he  scaled  in  like  manner, 
w 

The  scale  factors  for  the  material  properties  of  both  the  structure  and  the 
soil  have  a  great  hearing  on  the  achievement  of  usahle  scale  factors  for  the 
other  quantities.  The  ability  to  find  a  variety  of  simulants  to  provide  the 
combinations  of  scale  factors  necessary  to  bring  the  remaining  factors  within 
usable  limits  will  Increase  the  contributions  that  can  be  made  by  the  shock  tube. 


5.0  MODEL  MATERIALS 


In  order  to  properly  implement  the  scaling  laws  developed  it  is  necessary  to 
have  a  variety  of  model  simulants  available  so  that  usable  scale  factors  may  be 
obtained.  As  an  impetus  to  the  development  of  a  list  of  simulants,  a  variety  of 
some  materials  have  been  investigated.  These  include  prototype  materials  and 
possible  simulants  for  both  structure  and  soil. 

The  choice  of  the  proper  model  simulant  will  depend  on  the  total  scaling 
requirements  for  the  particular  experiment. 

A  rather  obvious  requirement  for  choosing  a  model  simulant  is  its  availability 
which  is  generally,  although  not  necessarily,  directly  related  to  a  reasonable 
cost  value.  Materials  not  easily  available  or  of  an  expensive  nature  would  not  be 
desirable  unless  their  use  is  unavoidable  and  the  results  justify  the  extra  effort 
and  expenditures.  Avoiding  these  choices  if  possible,  the  material  must  be  chosen 
whose  parameters  will  produce  usable  scale  factors  for  the  particular  case. 

One  important  fact  to  also  keep  in  mind  when  choosing  simulant  materials  which 
have  passed  the  criteria  mentioned  above,  is  the  ratios  between  model  and  prototype 
yield  or  ultimate  stresses.  These  ratios  must  not  be  less  than  the  scale  factor 
for  stress  being  used.  Whether  either  yield  or  ultimate  stress  values  are  employed 
in  the  choice  will  be  determined  by  the  particular  phenomena  being  investigated. 

For  example  it  may  be  desired  to  predict  stress  values  in  the  prototype  due  to  the 
shock  wave  impact  and  a  scale  factor  for  stress  of  1:10  is  indicated  in  projecting 
from  model  to  prototype.  If  the  prototype  is  constructed  with  10,000  psi  ultimate 
strength  concrete  the  model  simulant  must  have  an  ultimate  strength  equal  to  or 
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exceeding  1000  psi  or  the  ratio  of  ultimate  strengths  between  model  and  prototype 
must  exceed  In  this  example  case  the  yield  strength  need  not  follow  this  rule 

since  deflections  are  not  being  measured  and  only  the  breaking  point  is  of  import¬ 
ance. 


As  a  guide  in  formulating  a  testing  program  based  on  the  theoretical  back¬ 
ground  presented  in  this  report,  the  following  tabulated  list  is  provided.  Some 
"free  field"  experiments  in  the  shock  tube  should  produce  the  properties  not 
available  from  regular  known  sources.  Such  properties  as  were  available  from  the 
literature  investigated  have  been  included.  Blank  spaces  indicate  no  information 
available. 
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6.0  RECOMMEMDATIONS 


In  conjunction  with  the  shock  tube  facility,  the  scaling  law  derived  is  a 
tool  that  can  be  used  to  fill  voids  in  existing  knowledge.  The  particxilar  field 
of  underground  structures  subjected  to  nuclear  weapons  effects  has  many  such  voids. 
No  particular  approach  can  be  used  to  fill  all  of  these  in  one  fell  swoop. 

The  scaling  law  developed  will  provide  insight  into  some  of  the  major  voids. 
Other  deficiencies  in  present  knowledge  can  also  be  overcome  by  experimental 
work  based  on  this  law.  Some  of  the  areas  which  are  now  subject  to  question 
or  conflicting  views  can  be  filled  in  with  simple  yes  or  no  answers  that  can  be 
obtained  by  shock  tube  experiments.  For  example,  the  importance  of  damping  in 
the  structure  itself  has  been  debated.  The  primary  question  here  concerns  whether 
dEimping  within  the  structure  is  important  or  not.  Secondarily,  if  it  is  important, 
how  important  is  it?  Using  the  scale  law  (imposing  those  conditions  which  bear  on 
structural  deflection  and  vibration)  an  experimental  program  can  be  devised  to 
answer  the  primary  question.  If  two  model  cylinders  which  are  the  same  except  for 
their  dandling  characteristics  are  tested  \inder  the  same  conditions,  the  results 
should  provide  the  information  sought.  In  order  to  change  the  damping  character¬ 
istics,  one  cylinder  could  be  coated  on  the  Inside  with  a  damping  material  (such 
as  Vibrodamp).  If  damping  does  have  an  Influence,  a  quantitative  estimate  of  the 
degree  of  this  effect  can  be  obtained  by  a  series  of,  tests  with  varying  amounts  of 
damping  in  the  cylinder. 

Another  area  of  question  involves  the  importance  of  the  ratio  of  soil  mass  to 
structure  mass  in  the  case  of  shallow  burial.  Here  again,  by  varying  the  simulated 
soil  and  structure  materials,  different  vsilues  of  the  ratio  can  be  obtained.  A 
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wide  variance  in  two  tests  can  establish  whether  the  ratio  is  important.  If  it 
is  important,  more  subtle  variations  can  be  made  to  ascertain  more  quantitative 
information  as  to  the  effect.  Then,  the  depth  of  burial  can  be  varied  to  determine 
its  effect  on  the  problem. 

The  extent  to  which  some  questions  can  be  answered  depends  on  the  degree  of 
success  achieved  in  finding  suitable  model  simulants.  Simulating  the  soil  is  the 
greatest  problem  here. 

One  of  the  important  problems  that  can  be  answered  if  a  suitable  set  of  soil 
simulants  can  be  found  is,  "What  is  the  effect  of  the  unloading  slope  of  the  soil 
stress-strain  curve?"  If  two  simulants  can  be  found  whose  loading  characteristics 
are  similar  (can  both  be  scaled  to  same  prototype)  but  whose  unloading  character¬ 
istics  vary  to  a  significant  degree,  a  series  of  experiments  can  be  devised  to 
determine  the  effect  in  question.  By  the  seime  token,  the  Importance  of  the  depen¬ 
dence  of  the  effect  of  the  response  wave  on  soil  recovery  capacity  can  also  be 
determined.  Also,  the  effect  of  dissipation  in  the  soil  can  be  checked  in  the 
same  manner  (if  percyclic  damping  is  the  criterion). 

slight  variations  of  some  of  the  scaled  dimensions,  the  significance  of 
other  effects  can  be  observed.  The  effect  of  the  relation  between  the  arrival 
times  of  the  and  Cg  wave  and  the  natural  period  of  the  structure  can  be 
ascertained  by  a  series  of  experiments  where  everything  is  held  constant  except 
the  depth  of  burial. 

While  not  Included  in  the  scaling  law  as  derived,  the  asserted  existence  of 
viscoelastic  properties  of  the  soil  is  the  basis  of  a  problem  area.  The  existence 
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of  a  soil  simulant  which  exhibited  viscoelastic  properties  that  varied  rather 
widely  over  fairly  short  temperature  ranges  could  provide  qualitative  information 
that  would  shed  light  on  these  areas.  Such  a  material  could  be  used  in  a  series 
of  experiments  run  at  different  temperatures  and  the  results  compared. 

There  are  of  course  many  other  tests  that  may  and  should  be  performed  in  the 
shock  tube.  The  results  of  these  tests  should  then  be  made  available  to  theoreti¬ 
cians  so  that  further  advances  may  be  made.  With  new  and  revised  theories,  a 
revision  of  the  scaling  law  and  a  repetition  of  the  process  will  lead  to  more  refined 
and  exact  theories. 
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A  METHOD  OF  SIMDLATING  SOILS 
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APPENDIX  A 


A  Method  of  Simulating  Soils 

As  was  pointed  out  in  the  body  of  this  report,  one  of  the  most  difficult 
aspects  of  a  shock  tube  experimental  program  will  lie  in  the  simulation  of  soils. 

The  complicating  factor  here  is  the  non-linearity  of  the  stress-strain  curve.  In 
conjunction  with  an  internally  sponsored  AMF  project,  a  report  was  prepared  on  the 
subject  of  test  methods  for  underground  shock  problems.  One  part  of  this  report 
dealt  with  a  method  of  simulating  soils.  This  portion  of  the  report  has  been 
extracted  and  edited  and  is  included  in  this  appendix  as  a  possible  approach  to 
the  problem  of  soil  simulation. 

Since  this  report  covered  a  number  of  subjects  and  approaches  to  the  problem 
of  testing,  it  cannot  be  applied  directly  to  the  scaling  relationships  given  in  the 
body  of  the  present  report.  The  soil  'study  de&lt  with  pres  sure -specific  volume  lelaticai- 
ships  rather  than  stress-strain  relationships.  However,  the  method  given  here  will 
provide  an  insight  into  the  general  problems  of  soil  simulation  and  one  approach  to 
their  solution.  The  report  was  prepared  for  AMF  by  Mr.  Vincent  J.  Cushing  of  the 
Engineering  Physics  Company. 

A.l  Introduction 

If  one  is  going  to  employ  less  than  full-scale  pressures  in  connection  with 
soils  it  is  evident  that: 

1.  A  technique  must  be  found  for  manufacturing  a  soil  simulant  which  will 
provide  realistic  P-v  behavior  at  scaled  pressures; 

2.  The  P-v  behavior  of  real  soils  must  be  reasonably  well  understood  in  a 
quantitative  theoretical  manner  in  order  that  one  might  develop  a  cogent 
process  for  manufacturing  the  simulant  indicated  in  1  above. 
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In  our  discussion  here,  then,  we  will  concentrate  on  such  non-linear  stress-strain 
soil  characteristics;  discuss  a  theoretical  model  which  describes  real  P-v  behavior; 
and  thus  exhibit  the  parameters  which  must  be  controlled  if  one  is  to  fabricate  a 
soil  simulant. 

A. 2  Full-Scale  Soil  Behavior 

If  a  continuously  increasing  pressure  is  applied  to  soil-like  materials,  a 
pres sure -volume  (i.e.,  specific  volume)  relationship  is  generally  that  which  is 
indicated  in  Fig.  A.l.  Starting  with  the  initial  specific  volvime  v^,  the  P-v  curve 
is  initially  concave  downward  as  indicated  in  region  A.  In  this  region  the  (local) 
sound  speed  is  a  decreasing  function  of  pressure,  and  gives  rise  to  the  well-known 
dispersion  of  pressure  waves  in  soils. 

As  the  pressure  approaches  a  so-called  crushing  pressure,  P^,  the  P-v  curve 

flattens  out,  becoming  nearly  horizontal,  and  there  is  a  rather  sudden  decrease  in 

specific  volume.  In  this  region,  indicated  as  region  B  in  Fig.  A.l  the  soil  rapidly 

loses  its  structural  integrity,  becoming  essentially  fully  compacted  near  a  specific 

volume,  v',  as  indicated  in  Fig.  A.l.  As  the  pressxire  is  further  increased  this 

now  fully  compacted  soil  has  a  P-v  curve  as  indicated  in  region  C;  here  the  soil 

(  1) 

behavior  is  described  by  a  Tait  or  bulk  modulus  type  of  equation  of  state. 

If  now  the  pressure  is  unloaded,  the  soil's  P-v  relationship  will  return  along 
the  s£gne  Tait  type  of  curve;  and  in  fact  will  follow  such  a  curve  all  the  way  down 
to  zero  pressure,  finally  ending  up  with  a  specific  volume  of  v  ,  as  indicated  in 
Fig.  A.l.  In  this  overall  process  the  soil  has  become  fully  compacted,  and  in 
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PrcBsure 
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process  of  compaction  has  absorbed  from  the  pressure  source  an  amount  of  specific 

energy  as  indicated  by  the  hatched  area  in  Fig.  A.l. 

(2) 

Terzaghl  and  Peck'  indicate  that  for  certain  sands  and  clay  the  aforementioned 

crushing  pressure,  P^,  is  in  the  neighborhood  of  60-8O  psi.  On  the  other  hand, 

(•3) 

hovever,  the  Lavrence  Radiation  Laboratory'  '  indicates  that,  for  Nevada  tuff,  the 
crushing  pressure,  P^,  is  in  the  neighborhood  of  20,000-30,000  psi.  It  seems 
clear,  therefore,  that  we  must  be  analytically  and  experimentally  prepared  for  very 
extreme  variations  in  soil  characteristics. 

A. 3  Empirical  P-v  Relationships 

In  a  material  such  as  soil,  the  fine  details  of  a  P-v  relationship  depend  on 

(4) 

a  manifold  of  factors  such  as  strain- rate,  thermoviscosity,  porosity,  moisture 

(5) 

content,  etc.  However,  at  this  early  stage  of  our  analysis,  as  reported  here, 
we  will  limit  our  considerations  to  porosity- -and  we  shall  see  later  that  the 
characteristic  shape  of  a  soil's  P-v  behavior  can  be  described  by  consideration  of 
this  parameter  alone.  Since  we  are  not  going  to  consider  the  manifold  of  other 
factora,  we  shall  accordingly  limit  our  discussion  now  of  empirical  P-v  relation¬ 
ships  to  those  obtained  by  qtiasl-statlc  compression  wherein  strain-rate  and 
thermoviscosity  can  be  neglected. 

Region  C  of  Fig.  A.l  is  the  simplest  to  describe.  It  is  generally  called  a 
Talt  or  bulk  modulus  type  of  equation  of  state,  and  has  been  employed  extensively 
in  the  study  of  underwater  explosions.  The  most  notable  characteristic  of  such  a 
Talt  form  is  that  the  P-v  relationship  is  concave  upward.  Such  a  curvature  obtains 
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if  the  isentropic  bulk  modulus,  b,  is  a  linear  function  of  pressure 


b  s  -  v(SP/dv)g 

=  A(1  +  P/P^)  (A.l) 

where  A  is,  strictly  speaking,  a  function  of  the  entropy  S,  but  for  practical 
purposes  may  be  considered  as  a  constant.  From  Eq.  (A.l)  we  see  that  the  bulk 
modulus  increases  with  pressure.  Such  a  bulk  modulus  which  increases  linearly 
with  pressure  is  known  to  describe  quite  well  the  behavior  of  liquids  and  (already 
compacted,  i.e.,  without  porosity)  solids  at  very  high  pressures. 


If  the  isentropic  bulk  modulus,  b,  is  a  constant  (i.e.,  does  not  vary  with 
pressure) ,  then  it  can  be  shown  that  for  all  pressures  we  have 


c  =  hP/bp 


=  Pv  =  P  V 
o  o 


(A. 2) 


where  p=  l/v,  and  (P^,  v^)  is  any  point  on  the  P-v  curve.  Equation  (A. 2)  indicates 
that  the  product  Pv  is  equal  to  a  constant;  and  also  that  the  sound  speed,  c,  is 
constant  for  all  pressures.  This,  of  course,  would  be  the  case  for  a  constant 
bulk  modulus  or  constant  elastic  modulus — and  if  soil  materials  truly  were  typified 
by  such  a  characteristic,  our  problem  of  scaling  would  be  relatively  simple. 


Equation  (A.l)  indicates  a  bulk  modulus  relationship  for  which  the  sound  speed 
Increases  with  increasing  pressure  and  the  P-v  curve  is  concave  upward.  Equation 
(a. 2),  corresponding  to  the  constant  value  of  bulk  modulus,  corresponds  to  the 
situation  where  the  sound  speed  is  constant.  If  now  we  desire  a  relationship  where 
the  P-v  relationship  is  concave  downward,  evidently  we  require  a  bulk  modulus 
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relationship  which  is  a  decreasing  function  of  pressure. 


Empirically,  the  overall  curve  of  Fig.  A.l  can  be  approximately  described  by 


the  following  kind  of  a  bulk  modulus  relationship 

1  +  P/P 

» ■  *  [  1 "  (pA^)  1 


where  A,  P^,  and  P^  are  empirical  constants. 
Clearly  for  P  small, 

b«  A  [1  +  (P/Pg)] 


(A. 3) 


If  Pg  is  a  negative  number  ovir  bulk  modulus  is  the  desired  decreasing  function 
of  pressure,  and  we  would  therefore  have  the  P-v  curve  concave  downward  as  indicated 
in  Fig.  A.l  for  low  pressures.  On  the  other  hand,  for  pressures  sufficiently 
large  Eq.  (A. 3)  clearly  approaches  the  relationship  of  Eq.  (A.l);  and  therefore 
for  high  pressures  the  empirical  relationship  of  Eq.  (A. 3)  will  simulate  the  Tait 
regime  as  Indicated  by  region  C  of  Fig.  A.l. 


Equation  (A.  3)  is  a  three -parameter  (A,  P^^,  P^)  fit  to  an  empirical  curve 
typified  by  Fig.  A.l.  Clearly  expressions  with  a  larger  number  of  parameters  could 
be  devised  if  one  requires  a  better  fit. 


A. 4  Compression  and  Decompression--Irreverslbllity 

We  indicated  earlier  that  in  soils  the  compression-decompression  cycle  is 
Irreversible  in  that  the  soil  is  compacted  during  the  compression  phase,  but  does 
not  "decompact"  during  the  decompression  phase;  rather  the  decompression  portion 
of  the  P-v  curve  tends  to  follow  down  an  extrapolated  portion  of  a  Tait  form 
typified  by  region  C . 
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Today,  interest  in  underground  shock  is  centered  in  the  3,000-10,000  psi 
region,  i.e.,  ground  shock  induced  by  air  pressures  in  the  3,000-10,000  psi  region. 
Now,  if  the  crushing  pressure  of  the  soil,  P^,  is  large  compared  to  our  region  of 
interest,  then  in  actuality  a  typical  ccmipress ion -decompress ion  cycle  is  as  i..dicated 
in  Fig.  A. 2.  The  material  begins  with  a  specific  volume  of  v^,  and  is  compressed 
through  the  concave-downward  portion  of  the  P-v  curve  until  it  reaches  the  maximum 
pressure,  P^.  During  this  compression  portion  of  the  cycle  the  material  has  been 
partially  compacted  to  the  specific  volume  v^.  If  now  pressure  is  relieved,  the 
decompression  phase  of  the  P-v  curve  will  be  a  Tait  form  which  passes  through  P^; 
when  the  pressure,  P,  has  returned  to  ambient  the  specific  volume  of  the  material 
does  not  of  covirse  retiirn  to  its  initial  specific  volume,  v^,  but  rather  returns 
to  its  partially  compacted  specific  volume,  v^. 

A. 5  Theory  of  Compressibility  of  Porous  Media 

As  perhaps  suggested  in  the  previous  portion  of  this  Appendix,  our  requirement 
is  to  synthesize  a  material  wherein  the  bulk  modulus  varies  with  pressure  in  a 
prescribed  manner.  And  furthermore  since  we  may  be  interested  in  scaling  down  our 
pressures  by  a  large  ratio  (lO-KX)  or  more),  we  would  like  this  variation  of  bulk 
modulus  with  pressure  to  be  effected  in  the  scaled-down  pressure  regime. 

For  the  purposes  of  our  analysis  we  will  assume  that  we  are  given  a  desired 
P-v  relationship  such  as  that  indicated  in  Fig.  A.l.  We  shall  now  develop  a 
theory  of  material  compressibility  which  closely  duplicates  real  behavior;  the 
significant  parameters  in  this  theory  will  then  suggest  techniques  for  simulating 
a  medlvim  with  real  P-v  characteristics.  (We  will  not  attempt  to  describe  in  detail 
how  one  could  fabricate  such  a  soil  simulant--for  indeed  that  would  be  an  extensive 
undertaking  and  could  only  be  expected  within  the  framework  of  a  much  more  extensive 
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Pressure 


Specific  Volume,  v 

Figure  A. 2  i^essure  cycle  where  maximum  pressure  is  less  than  soil's  crushing 
pressure. 
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endeavor — but  we  will  make  some  suggestions  as  to  the  general  approach  which  might 
be  followed. ) 

A. 5.1  General  Picture — The  physical  picture  we  have  of  a  compactable  porous 
material  is  quite  similar  to  that  of  a  foam.  Generally  a  foam  is  made  up  of  a 
liquid  phase  with  an  inter-mixed  gaseous  phase,  and  perhaps  the  whole  system  stab¬ 
ilized  by  suitable  emulsifying  agents.  However,  for  a  solid  compactable  material 
we  think  of  the  foam  made  up  of  a  solid  material  which  is  interspersed  with  voids. 

In  real  soil-like  materials  these  voids  are  most  often  filled  up  with  air  and/or 
water;  however,  to  keep  out  picture  simple  at  this  stage  of  the  theoretical  develop¬ 
ment  we  will  assume  that  no  matter  is  in  the  voids,  and  hence  we  will  neglect  any 
compressibility  effects  of  the  matter  which  may  be  in  the  voids.  Later,  if  desired, 
the  theoiy  presented  here  could  be  embellished  in  a  straightforward  manner  to  account 
for  the  effect  of  matter  in  the  voids. 

If  the  system  of  voids  were  perfectly  regiilar  (i.e,,  if  the  voids  were  located 
with  something  like  crystalline  regularity,  and  each  void  had  the  same  volume  and 
shape),  we  would  have,  as  the  writer  sees  it,  a  very  imreallstic  sltiiatlon — and  the 
apparent  mechanical  properties  of  such  a  system  would  be  quite  different  (as  we 
shall  see  in  the  forthcoming  analysis)  from  the  P-v  behavior  of  observed  soil-like 
materials.  Rather,  we  believe  a  better  picture  is  one  wherein  the  voids  are 
randomly  distributed  with  respect  to  size,  and  we  assume  that  this  distribution  is 
imlform  throughout  the  volume  of  the  "solid  foam"--hence  for  analytical  purposes 
we  assume  that  the  distribution  function  describing  these  voids  is  invariant  from 
point  to  point  throughout  the  foam.  And  for  simplicity  in  analysis  we  will  assume 
that  this  distribution  is  Gaussian  with  the  mean  value  of  the  porosity  equal  to  a, 
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and  the  standard  deviation  in  porosity  equal  to  a  constant,  O  .  Let  us  now  look 

£l 

into  the  apparent  consequences  of  this  assvmed  physical  picture. 

A. 5 -2  Theory  of  Compaction — Let  the  specific  volume  of  our  foam  or  hulk 

material  be  v^,  and  let  the  specific  volume  of  the  solid  component  or  constituent 

be  equal  to  v  J  then  the  porosity,  a,  is  generally  defined  as 

a  =  1  -  V  /v,  (a. 5) 

s'  b 

Althoxigh  the  porosity,  a,  is  a  most  often  used  parameter  in  discussions  of 

compactable  material,  it  will  be  easier  for  our  analysis  to  employ  what  we  shall 

call  the  pore  ratio,  R,  which  we  define  as 

R  =  V  /v  (A. 6) 

p'  s 

where  v  is  the  specific  volume  of  the  solid  component  of  the  foam; 

8 

Vp  is  the  volvane  of  pores  per  unit  mass  of  the  solid  component. 

Now  the  specific  volume  of  the  bulk  material  is  equal  to 

V.  =  V  +  V  (a.?) 

b  8  p 

so  that,  using  Eqs.  (A. 5)  and  (A. 6),  the  relationship  between  our  pore  ratio,  R, 
and  the  normally  employed  porosity,  a,  is  given  by 

R  =  a/(l  -  a)  (A.8) 

Let  us  now  consider  a  sphere  of  our  foam  of  unit  mass  immersed  in  a  fluid 
so  that  a  hydrostatic  pressure,  P,  can  be  applied.  If  we  take  a  section  through 
a  diameter  of  our  sphere,  the  eu:ea  of  this  section  is  given  by 

A^  =7r^/^(3v^/4)^/3  (A.  9) 

and,  under  the  hydrostatic  pressure,  P,  the  total  force  acting  across  this  section 
is  given  by 

F  =7r^/3(3v^A)^/3p  (A.  10) 
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Now,  the  force,  F,  indicated  in  Eq.  (A.IO)  must  be  borne  by  the  load-bearing 
solid  component  of  our  foam.  Because  of  the  foam's  porosity  the  effective  cross 
section  of  load-bearing  material  is,  of  course,  less  than  the  of  Eq.  (A.9)j 
and  in  fact  the  effective  area,  A  ,  of  load-bearing  material  is  given  by 

B 

Ag  =7r^/^(3VgA)^/3  (A.ll) 

Let  P  be  the  stress  level  in  the  load-bearing,  solid  component  of  our  foam.  Then, 
s 

if  equilibrium  maintains  we  must  have 

8  8 

where  P  is  the  force,  as  indicated  in  Eq.  (A.lO)  which  is  acting  across  the  dia¬ 
metric  section  of  our  sphere  of  foam. 

Before  going  fairther  we  would  like  to  make  one  embellishment  to  the  relation¬ 
ships  that  have  been  derived  so  far.  We  should  like  to  assume  that,  at  each  point 
in  the  spherical  volume  which  we  are  considering,  the  solid  component  is  Indeed 

I 

load-bearing  as  long  as  the  stress  level,  P  ,  is  less  than  the  ultimate  stress 

0 

level,  P^,  of  the  solid  material  of  our  foam;  and  we  shall  assume  that  the  material 

is  no  longer  capable  of  carrying  any  load  whenever  the  local  stress,  P  ,  exceeds 

0 

P  .  At  the  Instant  that  P  exceeds  P  we  shsill  say  that  the  materleil  has  failed 
at  that  point;  and  we  shall  assume  that  the  stresses  in  the  solid  material  are 
continuously  redistributed  so  that  the  load  is  always  borne  by  the  material  which 
has  not  yet  been  crushed.  Taking  this  into  consideration,  then,  Eq.  (A. 12)  would 
be  modified  to  read 

p  -7r^/^[3v^(l  -  f)A]^'^^B  (A.13) 

where  f  Is  the  fraction  of  the  solid  component  of  our  foam  which  has  been  crushed. 
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Furthermore,  we  shall  assume  that  the  porosity  disappears  at  those  points 
where  the  load-bearing  material  has  failed;  that  is  to  say,  when  the  load-bearing 
material  fails  locally,  it  collapses  and  fills  the  voids  in  the  locale.  Taking 
this  into  consideration,  then,  Eq.  (A. 7)  and  Eq.  (A  lO)  must  be  modified  to  read 


V,  =  v  +  v(l-f) 
b  s  p'  ’ 


=  Vg  [l  +  R(l  -  f)] 


(A.l4) 


and 


[3Vg  [1  +  R(1  -  f))  A 


From  Eqs.  (A. 13)  and  (A. 15)  we  obtain 
Pg  =  p/g^/^ 


(A.15) 

(A.16) 


where 

G  =  [l  -  f  ]  /  [1  +  R(1  -  f)  ]  (A.17)  , 

Equations  (A.16)  and  (A. 17)  show  that,  at  (each  neighborhood  of)  each  point, 

the  stress  level,  P  ,  in  the  load-bearing  component  of  our  foam  is  a  function  of 
s 

the  local  values  of: 

a.  the  pore  ratio,  R; 

b.  the  failure  fraction,  f. 

We  shall  not  attempt  to  specify  explicitly  the  value  of  R  throughout  our  foam; 
rather  we  shall  be  satisfied  to  make  some  assumptions  about  the  statistical  dis¬ 
tribution  of  pore  ratio,  R.  First  of  all,  we  shall  assume  that  in  our  sphere  of 
foam  the  mean  value  of  the  pore  ratio  is  equal  to  R.  Accordingly,  the  average  value 
of  the  load-bearing  material's  stress  level  is  given  by 

(A.  18) 

where 

G  =  [1  -  f  ]  /  [1  +  R(1  -  f)  ]  (A. 19) 
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Next  we  shall  assume  that  the  pore  ratio,  R,  has  a  Gaussian  distribution  in 
respect  to  magnitude  and  that  the  standard  deviation,  a„,  is  a  constant. 

Since  R  and  hence  G  and  P  are  distributed  variables,  we  should  like  to  find 

s  ’ 

a  relationship  between  O^and  P  's  standard  deviation, cL  .  Before  we  can  find  such 

r\  S 

6 

a  relationship  we  must  first  express  P^  as  a  linear  function  of  R.  To  do  this  we 

shall  assume  that  throughout  our  foam  R  does  not  differ  too  much  from  Rj  and 

therefore  that  P  does  not  differ  too  much  from  P  •  Under  this  assumption  we  can 
8  8 

lineeurize  Eqs.  (a.16)  and  (A.I7),  i.e.,  simply  form  differentials  so  as  to  express 
AP  as  a  linear  function  of  AR, 

AP  =  (Pp/aK^A  R  (A. 20) 

To  a  first  order  approximation  Eq.  (A.20)  expresses  variations  in  P  as  a  linear 
function  of  the  variations  in  R;  hence  ^  6 ) 

Op  =  (2P/3)(G^'^  (A.21) 

The  distribution  function  for  a  Gaussian  distribution  with  mean  value 
described  by  Eq.  (A.18)  and  standard  deviation  described  by  Eq.  (A.2l)  is 
dP 

dX  =  exp-  [(P^  -  p3)V2ap  ]  (A. 22) 

Op  ® 

s 

Equation  (A. 22)  is,  of  course,  normalized  so  that 
00 

/  dX  =  1  (A. 23) 

—  00 

Using  Eq.  (A.22)  we  are  now  able  to  set  up  an  expression  which  ascertains 
the  fraction,  f,  of  the  load-bearing  material  which  has  failed;  that  is,  by 
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hypothesis,  the  fraction  of  material  wherein  the  local  stress  level,  P  ,  exceeds 

the  load-bearing  material's  yield  stress,  P^.  The  expression  for  f  is 

f=7p  dX  (A.24) 

u 

Equation  (A. 24)  can  be  expressed  in  terms  of  readily  available  tabulated  functions  as 
f  =  [  1  -  erf(xj]  /2  (A. 25) 

where 

2 

Xu  =  [(PyP)G  -  1  ]  /  [2/2  GO  p  ]  (A.26) 

(  7  ) 

where  we  have  made  use  of  Eq.  (A.21) ,  and  where  erf(x)  is  the  error  function  '  '  '  or 
error  integral.^  ^  Equations  (A. 25)  and  (A. 26),  then,  (together,  of  course,  with 
Eq.  (a. 17))  express  the  failure  fraction,  f,  as  a  function  of  the  applied  hydrostatic 
pressure,  P — and  in  terms  of  the  three  parameters,  P^,  R,  and  or^^.  Control  of  these 
three  parameters  will  provide  a  control  over  the  manner  in  which  failure  takes  place 
in  the  load-bearing  component  of  our  foam.  Figure  A. 3  is  a  plot  of  f  versus  P/P^ 
for  the  following  two  sets  of  conditions: 

a.  R  =  1,  and  Op  =  0.5; 

b.  R  =  0.5,  and  =  0.25. 

Once  the  failure  fraction,  f,  is  known  as  a  fvinction  of  the  hydrostatic 
pressure,  P,  we  can  utilize  Eq.  (A. 14)  to  ascertain  the  specific  volume  of  our  foam 
as  a  function  of  the  hydrostatic  pressure,  P — that  is  to  say,  we  have  now  determined 
the  P-v  relationship  for  our  theoretical  model  of  a  porous,  compactable  material. 

In  Fig.  A. 4  we  make  use  of  the  data  in  Fig.  A. 3  and  plot  P/P^  as  a  function  of  v^. 

The  solid  curves  in  Fig.  A. 4  indicate  the  behavior  of  our  compactable  material  if  the 
stresses  in  the  material  are  always  in  equilibrivan  with  the  hydrostatic  pressure. 

In  practice,  of  course,  when  such  a  material  is  subjected  to  a  hydrodynamic  shock 
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and  material's  yield  stress,  P^;  and  with  average  pore  ratio,  R,  and  standard 
deviation  in  pore  ratio,  cr  as  parameters. 


wave,  the  compactable  material  compresses  in  the  manner  of  the  solid  portion  of 
Fig.  A-4  until  the  collapsed  pressure,  P^,  is  reached — and  at  that  instant  the 
material  suddenly  collapses  and  fxilly  compacts  right  down  to  the  specific  volume 
of  the  solid  constituent.  Hence,  \inder  hydrodynamic  shock  the  observed  P-v  curve 
follows  the  dashed  portion  as  indicated  in  Fig.  A. 4^. 

A. 5. 3  Conclusion — In  subsection  A. 5.2  we  have  outlined  a  theory  of  compaction 
for  soil-like  materials.  Within  the  framework  of  this  elemej  cary  theory  we  have 
seen  that  the  significant  parameters  are: 

a.  P^,  the  ultimate  allowed  stress  in  the  solid,  load-bearing  component  of 
the  porous  material; 

b.  the  average  pore  ratio,  R,  throioghout  the  compactable  material; 

c.  the  standard  deviation,  a_,  in  the  distribution  of  pores  throughout  the 
compactable  material.  These  three  parameters  control  the  P-v  behavior  of  our  soil¬ 
like  material,  and,  in  particular,  control  the  hydrostatic  pressure  level,  P^,  at 
which  crushing  takes  place. 

In  this  first  approach  to  a  theoretical  model  we  have  assximed  a  Gaussian 
distribution  of  pore  ratio;  however,  further  consideration  would  perhaps  indicate 
whether  other  statistical  distributions  might  be  better  descriptive  of  the  real 
situation. 

But  whatever  distribution  function  is  employed,  the  fundamental  physical 
picture  which  we  portray  here  is  that,  with  increasing  hydrostatic  pressure,  P, 
there  is  a  probabilistic  progression  of  local  failure  in  the  load-bearing  component 
of  the  soil-like  material;  and  this  progression  of  failure  takes  out  of  play  a 
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fraction,  f,  of  this  load-hearing  material.  Consequently,  as  the  hydrostatic 

pressure,  P,  is  increased  the  effective  cross  sectional  area  of  the  load-bearing 

material  is  progressively  decreased.  As  this  process  is  continued  there  will  be 

reached  a  state  where  instability  sets  in;  the  fraction,  f,  rises  rapidly,  and 

the  soil-like  material  loses  its  Identity  and  collapses  (i.e.,  fully  compacts)  to 

the  specific  volume,  v  ,  of  the  solid  component  of  the  soil-like  material.  Once 

s 

the  porous  material  has  collapsed  or  compacted,  then  further  increase  in  hydrostatic 
pressure  will  find  a  P-v  variation  corresponding  to  the  bulk  modulus  equation  of 
state  which  is  proper  to  the  solid  constituent. 
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